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Abstract: We derive an integral expression for the leading-order type I-I-I three-point 
functions in the su(2)-sector of J\T = 4 super Yang-Mills theory, for which no determinant 
formula is known. To this end, we first map the problem to the partition function of 
the six vertex model with a hexagonal boundary. The advantage of the six-vertex model 
expression is that it reveals an extra symmetry of the problem, which is the invariance under 
90° rotation. On the spin-chain side, this corresponds to the exchange of the quantum space 
and the auxiliary space and is reminiscent of the mirror transformation employed in the 
worldsheet S-matrix approaches. After the rotation, we then apply Sklyanin’s separation 
of variables (SoV) and obtain a multiple-integral expression of the three-point function. 
The resulting integrand is expressed in terms of the so-called Baxter polynomials, which 
is closely related to the quantum spectral curve approach. Along the way, we also derive 
several new results about the SoV, such as the explicit construction of the basis with twisted 
boundary conditions and the overlap between the orginal SoV state and the SoV states on 
the subchains. 
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1 Introduction 

String theory was originally discovered as a natural field-theoretical formulation of the dual 
resonance phenomenological models of the strong interaction. Although this line of research 
was once abandoned after the advent of quantum chromodynamics, its basic philosophy 
is realized in a slightly different guise in the modern approach to the gauge and string 
theories, namely the AdS/CFT correspondence [1]. 

By now, a heap of evidence in support of the correspondence has been accumulated. 
Nevertheless, fundamental questions, such as how strings in AdS emerge as gauge-theory 
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collective excitations, are still left unanswered. To address such questions, it would be desir¬ 
able to establish non-perturbative approaches to analyze gauge theories. The integrability- 
based method, which is the subject of this paper, is one of such promising approaches. 

The integrable structure in the context of the AdS/CFT correspondence was first 
discovered in the spectral problem of planar AT = 4 super Yang-Mills theory (SYM) [2]. 
The subsequent rapid progress [3] culminated in the elegant non-perturbative formalism, 
known as the quantum spectral curve [4], which allows one to compute the spectrum at 
astonishingly high-loop order. Meanwhile, the integrability-based methods were extended 
also to other observables, such as Wilson loops [5-8] and scattering amplitudes [9-12], 
Lately much effort has been devoted to the study of three-point functions and structure 
constants [13-42], and a non-perturbative framework, called the hexagon vertex, was put 
forward quite recently [43]. Although powerful and remarkable, these non-perturbative 
frameworks rely on certain assumptions which have yet to be validated in gauge theories. 
The most notable among them are the so-called crossing and mirror transformations [44, 
45] . These transformations have their origin in the string world-sheet theory and are hardly 
visible on the gauge-theory side. Hence, to deepen our understanding of the duality, it is 
important to study the gauge theory more in detail and understand why and how such 
“stringy” characteristics can be borne out by the gauge theory. 

With such a far-reaching goal in mind, in this paper we revisit the computation of the 
leading-order three-point function in the so-called su(2) sector of AT = 4 SYM. A special 
class of such three-point functions (called type I-I-II or mixed in [41]) are well-studied in 
the literature and are known to be given by the scalar product between on-shell and off- 
shell Bethe states, which have a simple determinant expression. It was confirmed in [43] 
that the hexagon vertex also reproduces the same expression. On the other hand, a more 
general class of su(2) three-point functions (called type I-I-I or unmixed in [41]) are much 
richer in structure and the result is expressed by the complicated sums over partitions with 
the summand given by a product of three determinants. For this latter class of three- 
point functions, the hexagon vertex appears to be less effective and so far no closed-form 
expression has been obtained by that approach 1 . 

The main objective of this paper is to derive a new integral expression for such an 
intricate three-point function, with the hope of shedding light on its structure and setting up 
the foundation for future development. The method we employ is the so-called Sklyanin’s 
separation of variables (SoV) [46], which was previously utilized to study the scalar products 
(and the form factors) [47-50]. As illustrated in [50], to apply the SoV method to the 
periodic su(2) chain, we first need to introduce the twisted the boundary condition and 
then remove the twist at the end of the computation. Although such manipulation can be 
carried out straightforwardly in the case of scalar products, the removal of the twists turns 

1 Although no closed-form expression was obtained, the equivalence with the usual weak-coupling result 
was checked extensively by the case by case analysis [43]. 
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Figure 1.1. The three-point function corresponds to the partition function of the six vertex model 
with a hexagonal boundary. Here 0 1 ' 1 ' 1 = 6 1 ' 12 - > U 0* 31 -* etc. are the inhomogeneity parameters and 
u„ are the rapidities (n = 1,2,3). Blue and black lines correspond to the quantum space and 
the auxiliary space respectively. Here we are depicting the figure as if it is embedded in three 
dimensions, assuming that the angles between the blue and the black lines are 90°. There is a 
conical defect in the bulk with an excess angle 7r, which is in accordance with the hexagon vertex 
picture of [43]. 


out to be quite subtle for three-point functions. In order to circumvent this difficulty, we 
exploit the well-known correspondence between quantum integrable spin chains in 1(+1) 
dimensions and classical integrable statistical models in 2 dimensions [22], In the case at 
hand, the relevant statistical model is the six-vertex model and the three-point function 
turns out to correspond to the partition function with domain wall boundary conditions 
(DWBC) along the the hexagonal boundary depicted in figure 1.1. 2 If all angles are assumed 
to be 90°, there is a negative curvature defect with excess angle n in the bulk. 

The advantage of the six-vertex expression is that it makes manifest an extra symmetry 
of the problem, which is the invariance of the partition function under a 90° rotation. In the 
original spin-chain formulation, this cannot be seen easily as it corresponds to the exchange 
of the quantum space and the auxiliary space. Intriguingly, this symmetry is reminiscent 
of the mirror transformation employed in the non-perturbative approaches and we thereby 

2 The DWBC has been first defined on a rectangle by Korepin [51]. More generally, one can define them 
for any boundary consisting of 2n segments. The lattice with such a boundary has a curvature defect with 
excess angle (n — 2)7t. The case n = 1 was recently considered by Betea, Wheeler and P. Zinn-Justin in 

[52]. 
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Figure 1.2. Hexagon as the result of cutting the three-string world-sheet (in the sense of [54]). 
The world sheet has two point-like negative-curvature defects with excess angle n. If we cut in the 
temporal direction at the position of one of the defects, we obtain the hexagon with one curvature 
defect with excess angle n. This corresponds to the hexagon on figure 1.1 after being “flattened” 
to fit in two dimensions. 


call it the mirror rotation in this paper. 

The hexagon depicted in figure 1.1 can be thought of as the result of cutting the 
three-string world sheet along the temporal direction as shown in figure 1.2. In fact, we 
have encountered this hexagon configuration already in [25] for the EGSV configuration; 
in this case the contribution of the piece of the lattice associated with the excess angle 
factorizes and can be amputated, see also [26]. The rest of the lattice was brought to a 
rectangular form by the freezing trick and then evaluated as a scalar product. A similar 
procedure is at the core of the bootstrap method of [43] , where the three-point function is 
cut into two hexagons. In our case it sufficient to cut into a single hexagon, the second one 
degenerates into a Y-shaped junction of three seams. In general, cutting pants resembles 
the well known relation between closed and open string amplitudes [53]. 

The mirror rotation exchanges also the twists of the boundary condition and the global 
su(2) transformations acting on the spin chains. Importantly, such global transformations 
are always present for non-vanishing three-point functions. Thus, if we first mirror-rotate 
and then apply the SoV method, there is no need to introduce fictitious twists which 
will eventually be removed; the twists after the rotation are provided by the su(2) global 
transformations which exist already from the beginning. This feature allows one to express 
the three-point functions in terms of the SoV basis and the final result is found to be 

C 123 ~ d/u(x (a6) ; u a U u b ) d/al y a ; u a ) x (x (ab) - 0 (ab) ) 

J {ab) 

r(i(u+-u 6 -)) r(i-r(u+-x( afc )))r(i + i(u fe --x(^)) ( L1 ) 

r (i(y a - Yb)) r (l - i{Ya - x( a6 ))) T (l + i(y b - x( ab ))) 

X T(zi,z 2 ,z 3 ) 

where the product in the integrand is over the ordered pairs (ab) 6 {(12), (23), (31)}, 
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fi/i(x; w) is (up to a normalisation) the Sklyanin’s measures for the su(2) spin chain, 


<M x ; w) = ]~[ 

Xj^yi 


dxj A(x) A(e 2?rx ) 

27ri (x — w+)(x — w~) ’ 


A(x) c = Y[(xj - x k ), 

j<k 


( 1 . 2 ) 


and the factor T(zi, 22 , 23 ), whose expression can be found in section 4, takes into account 
the polarizations of the three states . 3 We also used the convention vr 1 = w ± i/2 as well 
as the shorthand notations of [18], namely a function of several sets of variables means the 
double product over all values of arguments, 


/(x)^ f n/(^ /( x -y)= n ( ls ) 

rrEx aiGx,yG y 

A notable feature of our result is that all the data characterizing the three operators, 
namely the rapidities u a and the inhomogeneities 6^ 12 \ 0( 23 ', $ 3 , appear only through 
the so-called Baxter polynomials , which in the convention (1.3) read 


Q 0 (««o(») = (® - ° {ab) ) , Qu a (x) = (x - Ua). (1-4) 


This feature would have two important potential implications. First, for a certain class 
of three-point functions, it is known that the one-loop result can be obtained from the 
tree-level result by judiciously making use of the inhomogeneities [21, 25]. Although such 
a method hasn’t been developed for a general class of three-point functions studied in this 
paper, our expression would provide an ideal starting point for such exploration since the 
dependence on the inhomogeneities takes a simple factorized form (1.4). Second, and more 
importantly, our result may provide some clues about how to utilize the quantum spectral 
curve approach [4] in the computation of the structure constants. The hexagon vertex 
approach [43], although non-perturbative, is only effective for sufficiently long operators. 
In order to study operators with finite size in full generality, it would be necessary to 
incorporate the method of the quantum spectral curve into the hexagon-vertex framework. 
Since the essential ingredient of the quantum spectral curve is the so-called P -p system, 
which is the finite-coupling analogue of the Baxter polynomials, expressing the three-point 
functions using the Baxter polynomial as in (1.1) may be regarded as a step toward such 
an ultimate goal. 

The rest of the paper is structured as follows. In section 2, the separation of variable 
for Heisenberg XXX 1 / 2 spin chain is discussed in detail. In particular, we derive explicit 
expressions for the SoV basis with twists at both ends of the chain, generalizing the result 
known in the literature [47-49]. In order to apply the SoV method to the three-point 
function, we then study how the SoV basis behaves when the spin chain is cut into two. 
We first derive a recursion relation obeyed by the overlap between the original SoV state 

3 A nice feature of this representation of the three-point function is that the homogeneous limit 6'' ab '> —>• 0 
is obvious and can be taken before performing the integral. 


- 5 - 




and the SoV states in the subchains, and then solve it utilizing the explicit expression for 
the basis. In section 3, we elucidate the symmetry of the domain wall partition function of 
the six-vertex model under the rotation by 90° and how it translates into the property for 
the scalar products of the spin chain. Of particular importance is that twists and the su(2) 
global transformations are exchanged under such a rotation. Then, in section 4, we derive 
a new integral expression for the three-point functions based on the results derived in the 
previous two sections. As briefly described above, the basic strategy is to first perform 
the mirror rotation and then apply the SoV method. We end with the conclusion and the 
future prospects. Several appendices are provided to explain technical details. 


2 Separation of variables for Heisenberg XXX 1/2 spin chain 

In this section, we construct the SoV basis for the XXX^ spin chain. According to 
Sklyanin’s recipe [46], the separated variables are the operator zeros of the B(u) operator, 
b(u) = b nt 1 (rt — Xk) • Together with the diagonal entries of the monodromy matrix 
A{xk) and D{xk) the separated variables xk can be used to construct sets of pairs of 
mutually conjugated variables. The separated variables are used as an alternative to the 
Algebraic Bethe Ansatz. We will denote by x = {xk]^ =l the eigenvalues of the separated 
variables {xfc}^ =1 and by |x) the corresponding eigenvectors, 

x fc |x) = Xfc|x) . (2.1) 

Since we can relate the SoV bases for chains with twists in different position using (2.18) 
and (2.19), it is sufficient to construct explicitly the basis for the left-twisted chain. The 
construction of the SoV basis for the anti-periodic chain [47] can be obtained as a particular 
case. The basics of the XXX x / 2 separated variables were described in [47, 50] and we refer 
to these works for more details. 

The main obstruction to construct the separated variables for the su(2) symmetric 
XXX 1 j 2 spin chain is that the B{u) operator is nilpotent and as such not diagonalizable. 
In order to apply the SoV formalism, one needs to introduce twisted boundary condi¬ 
tions, which breaks the su(2) symmetry in a minimal way and renders the B(u) operators 
diagonalizable. 


2.1 Twists 

The most general off-diagonal twist can be realized with an 5 1(2) matrix 


K=^° ^=e*“ aCTa , det K = ad — be = 1 , (2. 

where a a are generically complex numbers (real, if we consider a su(2) twist) and a a are 
the Pauli matrices in the auxiliary space. The twisted monodromy matrix Tk(u) is defined 
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Figure 2.3. Yang-Baxter 
equation for the twist matrix. 


Figure 2.4. s((2) invariance of the Lax matrix 
L (u), with K = e lrXaa the twist, i.e. rotation 
in auxiliary space, and gK = e 2z “ a s the corre¬ 
sponding rotation in the quantum space. 


by 

T k (» = KLi(u)---L l (u) = ( Ak ^ (2.3) 

K> \C K (u) D K (u)J f 

and it obeys the Yang-Baxter relation due to the su(2) invariance property of the R matrix 

Roo' (u) Kq K 0 ' = K 0 ' Kq Roo' ( u ), (2.4) 


with the index in Ko, Ko' representing the space in which the matrix K acts, as illus¬ 
trated in figure 2.3. This helps to show that the twisted monodromy matrix Tk(m) obeys 
the same Yang-Baxter equation as the untwisted matrix, and therefore its matrix ele¬ 
ments An, i?K) Ck, -Dr obey the same commutation relations as the non-twisted ones 
A, B, C, D. 

The property (2.4), which can be understood as su(2) invariance of the R matrix, is 
inherited by the Lax matrix L n (u) 


{<+aa) = e i aa K+a“) ^ 


(2.5) 


with the corresponding Pauli matrix at the site n of the spin chain. The Lax matrix 
for the XXXx /2 spin chain is given by 


Ln (u) 


(u + iS z n iS n A 

\ lS n U-iS^J ’ 


where S'" are the su(2) generators at site n, 


S a = \(<j a 1+ a a 2 + --- + ol). 


( 2 . 6 ) 


(2.7) 


The property (2.5) can be represented graphically as in the figure 2.4. 

The invariance property of the Lax matrix is also that of the untwisted monodromy 
matrix T(it) = Li (u )... L l(u), 

K T(u) K _1 = g^ 1 T(u) g K , g K = e 2iota sa . (2.8) 
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Introducing a twist has several consequences, notably changing the spectrum of the 
conserved quantities and modifying the expression of the B operator. The changes on the 
twist and on the monodromy matrix are correlated as follows 

K—>UK U -1 , T k (u) —»• UTk(u) U _1 (2.9) 

with T^(it) = gu T k (u) g(j 1 • 


The rotation in the auxiliary space, U, mixes up the elements A, B, C, D of the monodromy 
matrix, while the rotation gu affects only the quantum space. Since the conserved quantities 
are generated by the trace of the monodromy matrix, the spectrum of the twisted chain 
depends only on the eigenvalues of the twist matrix ( e lK , e _m ) via the twisted Bethe Ansatz 
equations 


L 


n 


Uj — Ofc ~b i/2 
Uj Ok i/2 


= e 2iK 


M 


n 


Uj — Uk + i 
Uj — Uk — i 


( 2 . 10 ) 


Let us now investigate the effect of changing the twist on the SoV basis. As we will 
show later, the left-twisted SoV basis will be constructed with the help of the raising-like 
operators 


Ak(u) = a A(u) + b C(u) , (2-11) 

and the basis will diagonalize the operators 

Bk(u) = a B(u) + b D{u) . (2-12) 


Any transformation of the twist which leaves the ratio a/b constant is therefore keeping 
the SoV basis unchanged. We conclude that the left SoV basis is left unchanged by the 
transformation 


K 




(2.13) 


The SoV bases are thus associated to the equivalence classes of twists under the transfor¬ 
mation (2.13). A representative for the equivalence classes can be chosen as 


K ~ 




or 


K ~ 


1 

\/l + 




= e Za g-| ln ( 1 + |z| 2 ) 0 - 3 e - 2 O-+ , 

(2.14) 


where (J± = (a 1 ± o' 2 )/2. The first choice has the advantage to be abelian under multipli¬ 
cation, while the second is unitary. 

The twists can be introduced at different positions of a spin chain. We can put the 
twist matrix at the left or right end, or at both ends of the spin chain 4 . The spin chains 

4 We can even put the twist in the bulk of the spin chain. 
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with these twists are called left-twisted, right-twisted, and double-twisted. We need to 
prepare the twisted chains to tailoring i.e. cutting the chains in two pieces each retaining 
a twist, so we will consider together the three types of twists. The twisted monodromy 
matrices are denoted as the following 


T Kl |o ( u ) =KiLiO)---L l (w), 
T 0|k>) =Li(u)---L l (u) K 2> 
T Kl|Ka («)=K 1 L 1 («)---L L ( U )K 2 . 


(2.15) 


where the twist matrix can be taken as any 2x2 complex matrix with unit determinant 

K 1= ( ai M, Ko=( a2 det K x = det K 2 = 1 . (2.16) 

\ ci d] J \c 2 d 2 J 

As the equation (2.8) suggests, the monodromy matrices with the twists in different po¬ 
sitions can be related to each other by rotations in the quantum space. For example, the 
right-twisted monodromy matrix can be written as 


To|k(«)=Li(u)---LlWK (2.17) 

= K (g K Li(tx) • • • L L (u) g” 1 ) = g K T K |o(«) g" 1 • 

It is then clear we can relate the SoV states for the left and right twisted spin chains as 
follows 


-l 

K 


l X )o|K Sk I X )k|0 ’ 0|k( X I K|o( X I§: 

and this relation can be generalized readily to the double twisted case 


=k 2 


/KiK 2 |0 


K 1 I k 2 K 1 K 2 |o( X I g; 


-1 

K 2 


(2.18) 


(2.19) 


2.2 Explicit construction of the SoV basis for the left twisted chains 

As explained in [50], the eignevalues of the separated variables, x, are related to the values 
of the impurities 0 = {0k)k=i by 

Xk — Ok T 2 Sfc i Sk — 1, k — 1,..., L . (2.20) 

For simplicity, we will denote alternatively the SoV basis by the values of the signs Sk, 

I x )k|o = |si,S2,...s l ) K | 0 , (2.21) 

with the obvious choice of signs according to (2.20). The number of + and — signs will be 
denoted by N x and N * respectively, with 

+ = ( 2 . 22 ) 

j =1 k =1 
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The right/left SoV basis can be constructed by applying sign-flipping operators to the 
the state with IV*= 0. We use that the diagonal matrix element J 4 K | 0 (u) of the twisted 
monodi'omy matrix with u = x k acts as a shift operator [50] 5 


■^k|o (xk ) \xi j • ■ • j x k > ■ ■ • 5 Tl)k|o Qo {x k ) \x\ x k x j/), 

K |0 (xi, ■ ■ • , Xk, ■ ■ ■, X L I A Kl0 (x k ) = (x k ) I K|0 (xi, ..., x k - i, . . . , XL I, 


with 


L 

Qe(x ) = - 0 k ), Q ± (x ) = Q(x±) = Q{x ± i/2). 

k =1 


(2.23) 

(2.24) 


(2.25) 


We can now construct the SoV ket-base and its dual bra-base starting from the reference 
states , —,... — )k|o an d k|o(+> ■ ■ • , +| respectively as follows, 


o 

II 

A K\o(^k ) 

_ Qo ) 

2 

y>, 

\i L ) = !-»■■■ ,-> K |0. 

K|0< X l = 

L 

n 

k =i 

i i 

§- 

1_1 

i ~ s k 

2 

1 

+ 

o 

a 

II 


(2.26) 


The identification of the reference states with | \. L ) and ("f L | can be done by noticing that 
they are eigenvectors of i? K | 0 (u) since 


-Bk|o( w ) = aB(u) + bD(u ) , 


and 


B(u) | l L ) = (t L \B(u) = 0 , D(u) | ) = Qt{u) || L ), (| L I D(u) = ( f L | Q~ e {u ). 

This, together with the relations (2.19) completes the construction of the SoV basis with 
two twists. 


2.3 Main results for the SoV basis 

Having an explicit realization of the SoV basis helps construct the main building blocks 
which are necessary to compute scalar products and correlation functions. In [50] some of 
these building blocks were determined from the functional (difference) equations they obey. 
But the difference equations do not completely fix the solution, so the initial condition had 
to be fixed in [50] by matching with some known cases. In this paper, equipped with the 
explicit construction of SoV basis, we are able to fix the ambiguities and determine all the 

5 Let us notice that the action of the operators ^4 k|o(^C fe) on the ket/bra SoV basis is the same as that of 
the right/left-ordered operators used in [50]. We are therefore going to skip the normal ordering sign. 
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relevant quantities for computing three-point functions. We present the main results with 
general twist and leave the derivation for the appendices. 

1. The measure. One of the basic property of the SoV basis is its completeness, and 
we will use often the resolution of identity 

1 = J2 l x ) Kl |K 2 M Ki ,k 2 ( x ) Kl |K 2 ( x l • ( 2 - 27 ) 

X 

The Sklyanin measure M Kl | K ,( x ) is nothing else than the inverse square norm of the or¬ 
thogonal SoV states 


k iIk 2 


(x'| 


X 


k iIk 2 


= /x 1 (x) <L 

^Ki|K 2 V > x 


S x ', x = S x[tXl ...S <tXL 


(2.28) 


In appendix A we show that it does not depend on the position of the twists and it is given 

by 


Mb 


x = 


n j<k( d j - e k)(Qj - Ok + i)(0j -0 k - i ) 


" k i | k 2 ^' (6i2 )L 

where by 612 we denoted the matrix element of 


Res [--- 

L(y - 0+ )(y - o ). 


K 12 = KiKa = 


012 &12 

C 12 d \2 


, det K 12 = 1 . 


. (2.29) 


(2.30) 


2. The vacuum projection. Another important ingredient for computing the scalar 
products is the projection of the SoV states on the pseudovacuum, /k 1 |k 2 ( x ) = Ki|K 2 ( x |t L 
). In appendix B we show that in the double twisted case, it is given by 


4,ik 2 ( x ) = . (2.31) 

The numbers N± are the numbers of pluses and minuses in the state k 1 |k 2 ( x I an d they are 
given explicitly in ( 2 . 22 ) in terms of the variables x. 

3. The splitting function. Let us consider a double twisted spin chain of length L. 
We cut the double twisted spin chain into one left twisted and one right twisted subchain, 
with length L\ and L 2) respectively, L\ + L 2 = L. The SoV basis for the double twisted 
spin chain can be related to the bases of the subchains as 


I x )k 1 ik 2 = 5 Z M Kl | 0 (yi) Mo,K 2 (y2) $(yi;y 2 |x) |yi) Kl ,o ® ly2) 0 , K2 - ( 2 - 32 ) 

yi>ya 

The overlap of the bases, $(yi; y 2 |x), or splitting function , is thus defined by 


^( yi ;y 2 |x) = Ki|0 (yi|® 0|K2 (y 2 | |x) Ki|Ka . (2.33) 

The splitting function obeys a set of difference equations for the variables on the subchains, 

&i 2 Qx(yi,j)^(y 1 ;y 2 |x) = h Q Y2 (yij) Qe 1 (yi,j)^(- • ■ ,yi,j-*,■■■ ; y 2 l x ) (2.34) 
bu Qx(y 2 j) d>(yi; y 2 |x) = b 2 Q y2 (y2j)Qg 2 (y2j)^{yi', ■ , m.j +*,••• |x). 


ll 






Here 0\ and 0 2 denote respectively the inhomogeneities of the left and the right subchains, 
e = 6 1 U 0 2 . These equations, derived in appendix C, have to be supplemented with an 
initial condition. Since the recurrence does not concern the variables x, we need a separate 
initial condition for each x. To this end, we consider a simple particular configuration and 
define 


^( x ) = Kl |o(+>-" +I® 0 |k 2 (-•••H x )k 1 |k 2 = (t Ll |®(t L2 |x) K i|K2 = (t L |x) Ki|K2 . (2.35) 
In appendix B we show that 

<Kx) = (~bi) N+ (b 2 ) N ~ ■ (2.36) 

The final result for < f?(yi;y 2 |x) after solving the difference equations (2.34) with the initial 
condition (2.36) is 


<f>(yi;y 2 |x) = twist x Gamma, 


twist = {ib 12 /b\) N + (-ibi 2 /b 2 ) N - (-h) N + (b 2 ) N ~ 


with N± = \L =f i^2(x k - 6 k ) N^ a = \L a =F i Y^(y a ,k - 9 a ,k 


k =1 


k= 1 


(2.37) 

(2.38) 

1 , 2 ), 


. r (get - e~ 2 )) r (i -<( yi - op) r(i + i( x - oj-)) r (1 -»(x-g 2 ~)) 

r(i(yi-y 2 )) r (l - z(y 2 - 0 2 )) T(l + z(x-yi)) T (1 — i(x — y 2 )) ’ 
where we used the shorthand notation (1.3). 

3 The scalar product as a six-vertex partition function: rectangle and 
rectangle in the mirror 

In this section we show that the scalar product between Bethe states has a different rep¬ 
resentation where the roles of rapidities and inhomogeneities are exchanged. We call this 
representation the mirror representation''. Upon this transformation, a global rotation is 
transformed into a twist and vice-versa. 

Our analysis is based on the six-vertex representation of the scalar product in terms of 
the Gaudin-Izergin-Korepin type determinant [51, 56, 58], found in [55]. This representa¬ 
tion has a remarkable symmetry under rotations with 90°. We can imagine the rectangular 
six-vertex lattice as a discrete world sheet of an open string, obtained by cutting the cylin¬ 
der along the time direction. Exchanging the magnon rapidities and the inhomogeneities 

6 It should be kept in mind that this in not exactly the same as the so-called mirror transformation in 
two-dimensional integrable held theories which transforms x it. 
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is like exchanging the space and time direction, hence the name we gave to this symmetry. 
Sewing the rectangle along the space direction, we obtain another with the space and time 
exchanged. 

After reminding the mapping between the pairing of Bethe states and the Gaudin- 
Izergin-Korepin determinant, we work out the correspondence between the transformations 
of the six-vertex configurations and the transformations of the monodromy matrix. The 
next step is to introduce rotations in the quantum space and twists in the auxiliary space 
and to show that they transform into one another under the mirror transformation. The 
two representations of the same object lead to two integral representations of the scalar 
product based on separated variables, one of them which appeared in [50] and the second 
being new. In the next section, the same techniques are used to write the three-point 
function in the mirror representation. 

3.1 The two-point function as (partial) domain wall partition function. 

The off-shell/on-shell scalar product of Bethe states can be written [55] in terms of a (par¬ 
tial) domain wall partition function, (p)DWPF [51, 57], and as such it has a representation 
in terms of the Gaudin-Izergin-Korepin determinant. A straightforward way see this rela¬ 
tion is to use the transformation property of the operators B{u) to C{u) which was proven 
in [42] 

B(u) = -E 2 C(u) S 2 1 , 

L 

with s 2 = i L J] crl = e inS2 

71=1 

The off-shell/on-shell scalar product is then 

(v|u) = (| L \C{v\)... C{y m )B{ui)... B(u m )\ | L ) (3.3) 

= (-ir<| L | B(vi )... B(v m ) S 2 B ( Ul )... B(u m )| t L > 

Using the Gauss decomposition in equation (3.6) below with £ = 7 t/ 2, as well as the highest 
height condition S' + |u) = 0 and the charge neutrality, one gets immediately 

(v|u) = (- l) m (i L \B ( V1 )... B{vm)B{u 1 )... B(u m )(S -) L - 2m I t L ) • (3.4) 

Let us now consider two Bethe states with global rotations, at least one of them, say 
the first one, being on-shell 


(3.1) 

(3.2) 


l^i) = g*K), * = 1,2. (3.5) 

A Bethe state rotated with an SU( 2) element can be labeled, as pointed out in [41], by an 
element of the coset space SU(2)/U(1). The coset structure appears because the rotations 
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with e iaS3 are acting trivially by multiplication with a phase. The generic element in 
SU(2)/U(1) can be parameterized as 

g z = e CS--CS+ = e zS- e -Mi+\*\')S* e -*s+ ) (3 .6) 

with z = C/ICI tan |^| and S± = S 1 ± iS 2 . When acting on an on-shell, highest weight 
Bethe state (5 + |u) = 0) the rotation can be brought close to the vacuum where it becomes 

g 2 |u> = (1 + \z\ 2 ) m ~ L / 2 e zS ~\u) = (1 + \z\ 2 ) m - L / 2 B( Ul ) . .. B(u m ) e zS ~ | ) (3.7) 

= {l + \z\ 2 ) m B(u 1 )...B(u m )g z \\ L ). 

Inspired from (3.3) we will represent the scalar products as pairings of ket states, (V2 11 i/^i) I ^2), 
using the su(2) singlet state, or two-vertex (V 2 I introduced in [41, 42], The explicit expres¬ 
sion of the two-vertex in the su(2) sector is 

<v 2 | = Y. (ei,...,£i|0 (£L,...,ei[S 2 . (3.8) 

ei,...,e i/ =t,4 

Using this formalism we get for the scalar product of two rotated states 
(V 2 ||V’i)|V’ 2 ) = (v 2 | g 2 |u 2 ) gi |ui) 

M 2 Ml 

=(-i) Ma a L 1 rt b m g 2 _i gi n 1 t L ) • (3-9) 

j =1 *=i 

The vertex (V 2 | helps to transfer spin chain operators from one chain into another following 
[41, 42]: 

(V 2 \A b \u) = (V 2 \D( 2 \u), (V2 \bW( U ) = -(V 2 \B ( - 2 \u), (3.10) 

(V 2 | DW(u) = (\ 2 \A^ 2 \u), (V 2 |C'( 1) («) = -(V 2 |C( 2 )(u). 

It can also be used to transfer the rotations from one chain to another 

(V 2 |g( 1 ) = (V 2 |(g- 1 )( 2 ). (3.11) 

This relation is a consequence of the singlet property of the vertex (V 2 |(S’^ 1 ^ a + S^ 2 ) 0 ) = 0. 
Without loss of generality we can set in equation (3.9) g = g^" 1 gi — g 2 

M 2 M\ 

<v 2 i iVh) 1^2) ^ a L 1 rt b m n b ko e * s ~ 1 t L ) (3T2) 

j= 1 2=1 

L-M M 2 Mi 

= (L-MV <i L in^d)II^K2)(g-) L - M |t L ) , 

1 >' j= 1 2=1 

where in the first line we have neglected a factor which can be reconstituted from equation 
(3.7). The last line is, up to the prefactor, the partial domain wall boundary condition 
partition function, pDWPF, with M = M\ + M 2 < L. The dependence on the global 
rotations is through the factor z L ~ AI . 
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U —*~ 


a(u, 9) b(u, 6) c(u, 9) 


Figure 3.5. The domain wall parti¬ 
tion function of six vertex model. No¬ 
tice that the configurations at the four 
boundaries are fixed. 


Figure 3.6. The six different ver¬ 
tex configuration, with their weights de¬ 
fined in equation (3.13). 


3.2 Direct and mirror representation for DWPF 

The domain wall partition function computes the partition function of the six vertex model 
on an L x L grid with domain wall boundary condition, as is shown in figure 3.5. 

The sum in the DWPF concernes all the configurations involving the six types of 
vertices shown in figure 3.6, the weight of each type of vertex being given by 

a(u, 9) = u — 9 + i/2 , b(u, 9) = u — 8 — i/ 2, c(u,6) = i. (3.13) 

The DWPF can be alternatively defined in the language of the algebraic Bethe Ansatz as 

L 

Z L (n\e) = (i L \llB(u k -,e)\t L ), (3.14) 

k =1 

where 9 = {9 1 , • • • , 9 £,} is the set of inhomogeneities and u = {ui, • • • , ul } are the magnon 
rapidities. The mapping from the six vertex configuration to the algebraic Bethe Ansatz 
language, explained at length in [22], is based on the interpretation of the six non-trivial 
vertex configurations (3.13) as the non-zero elements of the Lax matrix L(tt). The black, 
horizontal lines correspond to copies of the auxiliary space, while the blue vertical lines 
correspond to the quantum space. Conventionally, the Lax matrix acts from NE to SW. The 
DWBC configuration with all the arrows pointing upwards on the upper edge corresponds 
to the vacuum | f ^), the horizontal lines with incoming arrows correspond to operators 
B(u ) and the horizontal lines with outgoing arrows correspond to operators C(u). The 
symmetry properties of the Lax matrix are inherited by the six-vertex configuration, 

L n(u) = L n{u) t,to = a 2 n a 2 L n (u) (ala 2 ) -1 , (3.15) 

where f, to mean transposition in the quantum and in the auxiliary spaces respectively. The 
simultaneous transposition in the two spaces, followed by the conjugation with a 2 in the 


- 15 - 


























two spaces, amounts to the rotation of the corresponding vertex by 180°. The conjugation 
is necessary to keep the orientation of the arrows unchanged. 

When applied to the untwisted monodromy matrix, the simultaneous transposition 
reverses the order of the sites on the chain, as one can see as well from the graphical 
interpretation, 


T(u) f ’ t0 = L L (u) t,to ... L !(«)*■*“ = L L (u)... U(u) = f(u) . (3.16) 

In components, this means for example that B t (u ) = C(u), where the bar denotes the 
reversal of the order of the sites in the chain. Taking the transpose in the quantum space 
and using the fact that the DWPF is symmetric in the variables 9i,... , 9l, so that the 
order of the sites is irrelevant, one gets 

L L 

z L (n\0) = <t L | n £*(«*; 0)1 i L > = (t L l n c(u k -,e)\i L ). (3.17) 

k= 1 k =1 

The second equality in (3.15) translates into the following equality for the monodromy 
matrix 


T (u) = S 2 a 2 T(u) (Sad 2 )" 1 , (3.18) 

with £ 2 defined in (3.2). Written in components, this gives the relation (3.1) between the 
B(u) and C(u) operators. Given that the action of £ 2 on the vacua is £ 2 | ) = (—1) L | ) 

and £ 2 I \. L ) = | ) the DWPF can take the alternative form 

L 

Z L (u\6) = (t L \HC(u k -,e)\l L ). (3.19) 

k= 1 

The equality of (3.14), (3.17) and (3.19) expresses the invariance of the DWPF under 
reversal of the arrows. If one views the action of the monodromy matrix as an evolution in 
a (discrete) time, the rotation with 180° clockwise can be viewed as a PT transformation, 
and the transformations properties above as a CPT invariance. It is instructive to go back 
and interpret equation (3.9) in terms of six-vertex configuration. The second line can be 
written in terms of a rectangular six-vertex configuration. Now we insert the resolution of 
identity between the two rotations and apply transposition and conjugation with £ 2 in the 
first block, 

M2 Mi 

(_l)i-Ar a < 4,^! n ^ ( «2j) gjigi U B Mk)\t L ) = (3.20) 

j= 1 k= 1 

M 2 Mi 

= (£1,--.,£l| S 2 g 2 YlB(u 2 ,j)\t L ){£l,...,£L\gl n5(u i i fc)|t L ). 

ei,...,£ L =t,t j= 1 k= 1 
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Figure 3.7. Relation between six vertex model configuration and the spin vertex for¬ 
malism. The blobs symbolize the action of the operator S 2 : an extra factor (—1) if 
£=t- 

In the six-vertex model picture, the first line corresponds to the ordinary domain wall parti¬ 
tion (with the global rotations inserted in the middle) whereas the second line corresponds 
to the configuration where the lower half is rotated by 180° (see figure 3.7). The action of 
the operator S 2 is symbolized by the blobs in the middle of the lines. The effect of such 
blob is a multiplication by a factor (—1) if e =t- 

Less obvious are the transformation properties under rotation of the six-vertex configu¬ 
ration with 90°. By the space-time analogy above we can consider this transformation as an 
exchange of space and time, which in the context of integrable field theories goes under the 
name of mirror transformation. We are using this term here, although there might be dif¬ 
ferences with other instances of mirror transformations. Under the mirror transformation 
we get another domain wall boundary condition configuration with the inhomogeneities 
and rapidities exchanged as in figure 3.5. The vertices in figure 3.6 transform as follows 

a(u,9 ) —> b(9,u ) = —a(u,0) , (3.21) 

b(u, 9) —> a(9, u ) = —b(u, 9) , 
c(u, 9) —> c(9, u) = c(u, 9). 

We conclude that the rotated DWPF shown in figure 3.8 is equal to the non-rotated one 
up to a sign. To compute this sign we need to know the parity of number of vertices of 
type a and 6, N a and Njj. Due to the particular type of boundary condition we consider, 
on each line there should be one type c vertex. For the DWPF the sign is then given by 

(-1 ) N «+N b = (_ 1 )L(L-1) = 1 (3.22) 

so that 

Z L (6 |u) = Z L { u|0) . (3.23) 
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Figure 3.8. We rotate the configuration in figure 3.5 by 90° clockwise, the 
expression for the partition function is invariant. The interpretations of the 
two diagrams are different. The roles of rapidities and inhomogeneities are 
exchanged in the two configurations. 

Therefore, the DWPF in the mirror representation is given by 

L L 

Zl(0\u) = (t L | n C(e k -u)\i L ) = U L | n B(0 k - u)|t L ). (3.24) 

k= 1 k =1 

3.3 Global rotation and twist 

The next step is to find the mirror representation of the partition functions in the six vertex 
model in presence of a global rotation. We suppose that, as in equation (3.12) the rotation 
acts directly on the vacuum. The partition function we study is depicted in figure (3.9). 
On the l.h.s. the partition function corresponds to 

lhs = ('t L \g t C{u 1 )---C(u M )\i L ) = (-l) L - M (i L \g- 1 B(n 1 )...B(n M )\t L ) 

= (l L \B( Ul )---B(u M )g\t L ), (3-25) 

where in the last two equalities we have used the conjugation with the matrix £2 and 
transposition in the quantum space respectively, as discussed in the previous subsection, 
and g _1 = £2 gt • After rotation with 90° the partition function acquires a sign 
(— l) Na+N b = (— l) M ( i_1 ), cf. equation (3.21), and the rotation in the quantum space 
g is replaced by the twist matrix K* acting in the auxiliary space 

rhs = (-1) M ( L - 1 )(| M |S Kt | O (0i) • • • t M > (3.26) 

We conclude that after the mirror transformation the rotation in the quantum space g is 
replaced with the twist K*. In particular ' 

g = e ,s- ™ K f = e za+ . (3.27) 

Since in the mirror representation the 13-operators are twisted, we can apply the SoV 
formalism. 

'In what follows, we will denote K* by K to avoid cumbersome notations. 
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Figure 3.9. The mirror representation in the presence of global rotation. After turning 
the diagram 90°, the global rotation plays the role of a twist. 


3.4 Two dual integral representations for the scalar product 

Both the expressions (3.25) and (3.26) can be written as multiple integrals using SoV 
method. For (3.25) the integral representation was derived in [50], after introducing an 
extra twist which can be subsequently set to zero. Since we are going to compute scalar 
products for sub-chains, the rapidities are not assumed to be on shell. Nevertheless one 
can apply the argument of [50] to carry on the computation. 

Consider first the lhs, eq. (3.25), which we denote using the same notation as above 
(but with |0| > |u|). In order to go to the SoV representation, we introduce a left twist 
K e |o and then take the limit e —> 0. For simplicity we take g = e zS . Only one term in the 
expansion of the exponent survives and the result is 

lhs ->• (i L \B e (ui)... B e (u M ) e zS |t L ), (3.28) 


In this way we represent this expectation value as a multiple integral in the separated 
variables. The derivation is the one from [50], after noticing that only one term in the 
expansion of the exponent survives (the one which compensates the extra S z charge L — M ). 
We give only the final result, 


lhs 


L-M 


(. L-M)\ 



dxj Q u (xj) e 2 ^- 1 )^ 
2ni Qq (xj)Qg (xj ) 


L 

Y[(xj-x k ) 

j<k 



(3.29) 


where the function Eg is defined by 


FI j<k(9j ~ °k)(0j - 0 k + i)(dj - 6 k -i ) 
TIj<fc (e 2n9 i — e 2?r6,fc ) 


(3.30) 


and the contour Cg encircles the sets G + and 6 . Now we write integral representation for 
the rhs, eqn. (3.26), 


rhs = (-1)^-1 ) iM Z L-M 



dxj 
2iri 


M M 

n ( x J-x k ) n 

j<k k =1 


Qe{x k ) e 2 < k ~^ Xk 
Qu {%k) Qu {xk) 


(3.31) 


where H u is defined similarly as in (3.30) and the contour C u encircles the all points u + 
and u“. 


19 














































4 Three-point functions in the SoV representation 

In this section, we compute the structure constant using the spin vertex formalism in 
[41, 42, 59]. After the mirror transformation, we can compute both the spin vertex and 
the wave functions of the external states in the SoV representation, and combining them 
we get the final result. We can of course do the computation with the generic twists, but 
in fact it is enough to consider the triangular twists of the following form 

K Za =e z “° + =(^ Z A , a = 1,2,3. (4.1) 

The reason is that the twists in the mirror representation come from the global rotations 
in the original representation. As is shown in (3.7), if we start with on-shell Bethe states 
which satisfy highest weight conditions, the most general global rotation can be reduced to 
the rotation of the form g z = e zS , together with some factors which are not relevant. So 
we can define our external states as \i/j) = e zS |u) without loss of generality. The global 
rotation e zS , upon performing the mirror transformation, turns into the twist of the form 
(4.1). 

4.1 Spin vertex and the mirror transformation 

We are interested in computing the three-point function for three operators belonging to 
the so-called left su( 2 )^ subsector of the so(4) ~ su( 2 ) L <g> su( 2 ) R sector of the AT = 4 
SYM theory. This sector is made by the scalar fields X, X, Z, Z which belong to the 
bi-fundamental representation of su( 2 )l <S> su( 2 ) r , 

\Z) = |t) L ® |t) R = Itt) , \Z) = ||> L ® ||> R = |U> , (4.2) 

l*> = lt) L ® |f) R = \U), \X) = -|f) L ® |t) R = -lit) ■ 

Again, we are using the vertex formalism from [41, 42, 59] to compute the overlap of the 
three spin chains. At tree level, the three vertex (V 3 I is composed by the three singlets 
(ij)(V 21 corresponding to the bridges connecting the piece ( ij ) of the chain i with the piece 
( ji) of the chain j. 


(V 3 I = (12)(V2|® (23) C^21 ® (31) {V 2 1 • (4-3) 

We consider the case where the non-trivial magnon excitations belong only to the left 
sector 8 . The three external Bethe states with global rotations g Za <g> g Za (a = 1,2,3) are 
given by 


®o) = Ha) h ® | $a) R 


(4.4) 


This class of three-point functions are called type I-I-I or unmixed in [41]. 



with the right components being rotated pseudovacua, 

l^a) L =gJu a ) L , |^a) L = gij t La ) R (a = 1,2,3). (4.5) 

The spin vertex also splits into two identical parts (V 3 I = L (V 3 | (g> R (V 3 |. This insures the 
complete factorization of the left and right sectors, 

<v 3 | |*i> 1*2) 1*3) = (<v 3 | 1^1) 1^2) |^ 3 » l ((V 3 |g, 1 |t Ll )gi 2 |t i 2 )gi 3 lt i3 )) R - ( 4 . 6 ) 

The right piece is easily calculted and is equal to ( 5 i 2 ) il 2 (^ 2 3 ) L 23 (^ 3 i) i31 , where L a b = 
\{L a + Lf, — L c ), (a, 6 , c = 1, 2, 3). 

As mentioned before, we take the following external states 

M a 

| il’a) = e ZaSa |u a ) = Yl B{u aJ ) e ZaSa I t L “), a = 1, 2, 3. (4.7) 

j =1 

We are going to concentrate from now on on the structure constant in the left sector and 
drop the L index on the states and on the vertex 

C , i L 23 = (V 3 ||V’i) l |V' 2 ) l |^ 3 > l . (4.8) 

The structure constant above has a representation in terms of 6 -vertex model partition 
function on the diagram of figure 4.10. Apart from the blobs, this diagram is nothing but 
a redrawing of the hexagon DWPF in figure 1.1 (in the hexagon DWPF there was no need 
to put blobs along the “bridges”). 


031 012 0 12 023 023 031 



U! U 2 U 2 U 3 «3 “I 



Figure 4.10. The three-point function before and after the mirror transformation. One can see 
that the role of impurities and rapidities are exchanged, as in the case of scalar products. In 
addition, the global rotation becomes twist in the mirror representation, as in the scalar product. 

We will perform a mirror transformation on the six-vertex configuration of the three- 
point function as is shown in figure 4.10. We rotate the left subchains clockwise and 
the right subchains anti-clockwise and combine together the neighboring subchains, as is 
shown in the shaded region of figure 4.10. Potentially, there is some minus sign coming from 
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the transition from the direct vertex (V 3 I to the mirror vertex (V 3 I. But since the signs 
coming from “bridges” in the direct and the mirror vertex are related, since ones constitute 
boundary conditions for the others, we neglect any overall sign which may occur. From 
the 6 -vertex configuration, we see that we need to compute the following mirror structure 
constant 9 


Q L 23 = <5i L 23 = <V 3 ||^ 12) >|^ 23 )> |^ (31) >. (4.9) 

In what follows, we put a tilde on the operators which are in the mirror representation. 
The three mirror external states are now 

\ft ah) ) = II B Zal _ Zb (0[ ab) ) ( ab ) = (12), (23), (31), (4.10) 

k=1 

where the operators B Za \_ Zb (9) are the matrix elements (first row and the second column) 
of the following double twisted monodromy matrices 

M a M b 

f Zal _ Zb (9) =K Za n L n (9 - u ad ) n L k (9 - u btk ) K_ Zb . (4.11) 

j=1 k=1 

We can construct the SoV states for the three double twisted spin chains, which are denoted 

as \*)z a \-z b - 

4.2 Spin vertex in SoV representation 

In order to compute the spin vertex, we apply the important property (3.10). In the 
presence of the twists, this property is modified to be 

<v 2 | <U(») = (V 2 | r>2U(»), = -(V 2 I (4.12) 

(V,|flSLW = (V 2 |<l-„(»). <V 2 |cfIU>) = -(V 2 I 

By putting one of the twists to zero, we obtain similar relations for the left or right twisted 
monodromy matrices. Notice that the left twisted monodromy matrix is translated to a 
right twisted one by the spin vertex and vice versa. Using these relations, one can show 
that 


(v 2 | |x)*|o|y)o|-* = do(y|x)z|o = (/U|o(x)) ^x,y (4.13) 

where <5 Xi y = 5 Xltyi ■ ■ ■ S XN , yN and ^~| 0 (x) is the Sklyanin measure. 

We can write the three-point spin vertex in the SoV using the resolution of identities 
of the SoV basis. Denoting by x.( ab ^ are the SoV variables associated with the state 

u Here “=” means equal up to some overall minus sign. 
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we write the 3-vertex as a triple sum 


(v 3 | = Y, (4.14) 

x(^) ,x(^) ,x(31) 

x z 1 |-- 2 (x (12) | <8> Z2 |_ 23 (x (23) | <g> 23 |_ 21 (x (31) |, 
where the coefficient function V(x( a6 )) is given by 

V^W 1 )) = (V 3 | |x( 12 )) 2l| _ 22 ® | x ( 23 )) 22| _ 23 ® | x ( 31 )) 23| _ 21 . (4.15) 


In order to compute this function, we first split the SoV states as is described in section 2.3 
appendix C, 

|x (12) ) 2l |_ 22 = X] ^do(yi) Mo|- 22 (y 2 ) , h(yi;y 2 |x (i2) ) |yi) 21 | 0 ® |y 2 > 0 |- 22 , (4.16) 

yi,y 2 

and similarly for |x (23 ^) 22 |_ 23 and |x^ 31 )) 23 |_ 21 . Acting the states on the three-point spin 
vertex and using (4.13), we obtain 

v(x (i2) ,x (23) ,x (31) ) = Y M*i|o(yi) ^ 2 |o(y 2 ) M«3io(y3) (4.17) 

yl.y2.y3 

X 4>(yi;y 2 |x (12) ) <h(y 2 ;y 3 |x (23) ) <h(y 3 ;yi|x (31) ). 


The eigenvalues for the three sets of separated variables of the subchains are 

yi = ui±§, y 2 = u 2 ±§, y 3 = u 3 ±f, (4.18) 

where y a = u a ± | is understood as y a ^ k = u ajk ± (k = !,-■■ , M a ). 


4.3 The wave functions 

The external states in the mirror representation are given in (4.10). Acting these states on 
the spin vertex (4.14) and using the property of the SoV basis 

Zi\—Z2 ( x l^2i|— Z2 (^) (Al — z %) Qx($) Zi \— 2 2 (x|; (4-19) 

we obtain the wave functions of the mirror Bethe states 

M a -\-Mb 

,.|-„A“‘V ( “* ) >=(-1 /«.!-„(x<“‘>) n %<.‘)(4“'’ ) ). ( 4 - 2 ») 

k=1 

where the phase factor (— l^ L o.b( M a+M b ) comes f rom the rewriting of products of Q-functions 

L a b M a -\-Mb 

n Qx a {of b) ) = (0 {ab) - x {ab) ) = (- 1 )^( M “+ M «) Y[ Q g(ab) (x^). (4.21) 

j=1 k=1 
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The functions / 2 i_ 2 i) (x) are the projection of SoV basis on the pseudovacuum and are 
given by 


=z„|- I »(x < “ 1 ' , lt"“ +M ‘). (4.22) 

From (2.31) the projection of SoV basis on the pseudovacuum for the triangular twists of 
the form (4.1) is simply 

2 a |_, 6 (x|t L ) = l. (4.23) 

4.4 The final result 


We can now assemble the results from the previous subsections and write down the final 
result for the structure constant. Plugging (4.14), (4.15), (4.17) and (4.20) into (4.9), we 
obtain 

c^=u^) LaMa U^ b E E 

a=l (ab) {x( ai, )}{y a } 

3 M a +M b 

x Il^io(y-) n^l-.( x(a6) ) $ (y^;yb|x (ab) ) Yl Qo^ b )( x k b) )i ( 4 - 24 ) 

a= 1 (ab) k= 1 


where we have defined z a b = z a — Zb■ In (4.24), the summation is over all the possible 
eigenvalues of all the SoV variables. The structure of the summand is the product of 
Sklyanin measures of the three spin chains and three subchains, the three splitting functions 
which originate from cutting the spin chains, and the three wave functions. 

Similar to the scalar product, there’s another way of writing C \‘ 25 which converts the 
sum over eigenvalues of separated variables to a multiple contour integral. The derivation 
of the integral representation is analogues to the scalar product in [50]. The result is given 
by the integral 


C 123 = factor X 


{ui,2,3±*/ 2 } 


n 

(a,b) 


r (»( u a ~ Uft )) r (1 ~ »( u a ~ x(afc) )) r (1 + ~ x(ab) )) 

r (*(y a - y b)) r (i - i( y a - x( a6 ))) r (1 + *( y b - x( a6 ))) 


x ( X (« 6 ) _ ffW) T(z h z 2 ,z 3 ) dvixW) d M (y a 


(4.25) 


where (ab) = ( 12 ), (23), (31). We denoted by T(z i, z 2 , 2 : 3 ) the factor in the integrand which 
depends on the twists: 


with 


T(Z!,Z2,Z 3 ) 


n 




(-Z a 




(4.26) 


Nl a 


Ma 

2 


M a 

"F ^ — u a,j)> 

3 = 1 


Nl = 


Ma + Mb 

2 


Ma+M b 

T l E ( X< f b) ~ U ^3 

3 = 1 


(4.27) 
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The measures cfyr(x( a6 )) and /r(y a ) for the spin chains and subchains are defined by 



(4.28) 



(4.29) 


Finally, the overall factor reads, up to a phase factor, 

1 


(4.30) 


factor = 


(Mi + M 2 )!(M 2 + M 3 )!(M 3 + Ml)! Mi!M 2 !M 3 ! 

x X\( z ab) Lab ~ Ma ~ Mb X\{z a y Ma 


L a b M a 



X 'U1UU2 ^ •U2UU3 X ^U 3 UUl X l — J Ul X '^U2 ^ —-U3 • 


The factor comes from the product of the common factor in (4.24) and the part of Sklyanin 
measures which do not depend on SoV variables, Eq. (3.30). 

As proved in [41], the z a dependence is completely determined as follows by the Ward 
identity when operators are primary: 

cf 23 oc ( Z12 ) L 12- M 1- M 2+M 3 ( Z23 ) l 23-m 2 -m 3 +m 1 ( z ^l 31 -m 3 -m 1 +m 2 _ (4. 31 ) 

One way to see this z-factorisation explicitly is to first compute the integral, use the Bethe 
equations and eliminate z’s 10 . Although we haven’t succeeded, it would be much more 
desirable if we could rewrite the integrand (using Bethe equations or Baxter equations) 
in such a way that the ^-independence becomes manifest. We leave this as an important 
future problem. 

5 Conclusion and prospects 

In this paper, we derived a new integral expression for three-point functions in the su(2) 
sector of AT = 4 SYM using Sklyanin’s separation of variables. In order to apply the SoV 
method, we first mapped the three-point function to the partition function of the six-vertex 
model with a hexagonal boundary and then performed 90° rotations, which we call mirror 
rotations. The SoV approach can be readily used after this manipulation without the need 
of introducing boundary twists. The intriguing feature of our result is that the rapidities 
(and the extra inhomogeneities) enter only through the Baxter polynomials, which are 
considered to be intimately related to the quantum spectral curve approach [4], In this 

10 Such manipulation was performed in appendix L of [43] in order to compare the predictions from the 
hexagon vertex with the weak coupling result. 
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sense, the result obtained in this paper may be regarded as a (small) step toward the 
quantum-spectral-curve approach to the structure constants. 

We focused only on the su(2) sector in this paper. It would be an interesting future 
problem to extend the analysis performed here to other sectors. Of particular interest 
would be the generalization to the sl(2) sector, where the SoV basis already exists [60]. 
In the st(2) spin chain, the quantum space and the auxiliary space belong to different 
representations in the conventional formulation. Thus, the mirror rotation, if it exists, will 
take a very different form. For the scalar products, it is possible to write down an integral 
expression 11 akin to the mirror representation in the su(2) sector [61]. It would be worth 
investigating if a similar expression can be obtained also for the three-point function itself. 

Another interesting future direction is to understand the relation between the mir¬ 
ror rotation employed in this paper and the “genuine” mirror transformation used in the 
worldsheet S-matrix approaches [44, 45]. At the moment, it is not clear (at least to us) 
whether such a connection exists at all. However, it would be very intriguing if we could 
establish the relation as it may pave the way toward understanding the string world-sheet 
theory from the perturbative gauge theory. Also for this purpose, studying other sectors 
will be useful 12 . 

In order to see if our expression has a neat finite-coupling analogue, it would be im¬ 
portant to study in detail how the expression would be modified at one loop. To this end, 
it would be helpful to apply the approaches developed in [21, 25], where it was shown that 
a certain class of the one-loop structure constants can be obtained from the leading order 
result by a clever use of the inhomogeneities. It would also be interesting to try to factor 
out the dependence on the angles z a b in (4.25) using the Bethe (or Baxter) equation and 
further simplify the final expression. 

In deriving the integral expression, we obtained several new results about the SoV 
basis, such as the explicit expression of the states in the presence of twists and the splitting 
function which determines the overlap between the SoV state in the original chain and the 
SoV states in the subchains. These results may be useful in other problems, such as the 
computation of the form factors [39, 62-64] and the entanglement entropy [65] in the spin 
chain. It would also be interesting if we could use such results to study the hexagon vertex 
in the SoV basis. This may give some hints about how to incorporate the quantum spectral 
curve techniques into the hexagon-vertex framework. 

We hope that the materials studied in this paper will play a foundational role in the 
future progress and help unravelling still enigmatic features of the gauge/string duality. 

11 The ordinary SoV integral expression for the scalar products in the sl(2) sector is already known in the 
literature [30, 60]. 

12 To address such a question, it might be helpful to formulate the mirror rotation used in this paper as 
some kind of the (anti-) automorphism of the underlying algebra, as was the case for the crossing and mirror 
transformations of the worldsheet S-matrix [44, 45] . 
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A The Sklyanin measure 


Using the explicit representations of the SoV basis derived in (2.26), let us determine the 
Sklyanin measure. Here and throughout the appendices we consider general sl(2) left and 
right twists with matrices Ki and K 2 with the following notations, 


K,; = 



KiK 2 = 



(A.l) 


aia 2 + biC2 a\b 2 + b\d 2 \ 
cia 2 + dic 2 ci& 2 + d\d 2 ) 

The Sklyanin measure is given by the inverse of the norm of the SoV state as follows 

( x '|x) = Mkj| K a ( x ) <5x',x • (A-2) 

In [50], the measure was shown to satisfy a certain difference equation and was determined 
up to the overall constant by solving the difference equation. However, using the repre¬ 
sentations of the SoV basis (2.19) and (2.26), one can derive the following more explicit 
formula 


^Ki|K 2 ($! + 2 S 1> • • • i + 2 Sl ) ~ (1 ill 


k L 


^KiK 2 |0 (8k ~ |Sfc) 

QoiPk ) 


II 1 


(A.3) 


By decomposing A KlK2 | 0 in terms of the A, B, C and D operators for the untwisted 
monodromy and using the conservation of the su(2) spin, the right hand side of (A.3) can 
be further simplified as follows 

C(9k - §«fe) 


Mki|k 2 (^! + | s i , • ■ •, + |sl) — ( bi 2 ) L (t L I 


k 




II 1 


(A.4) 


Note that the right hand side of (A.4) is nothing but the domain wall partition function 13 
Since the measure factor is known to satisfy the difference equation, in order to determine 


13 Indeed, using the su(2) symmetry, one can write the right hand side of (A.4) alternatively as 

B{6 k - fsG 


(biA L (l L 


n 


Qe( 0 k — isk) 


ItY 


(A.5) 
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it unambiguously, it is enough to calculate it at one particular value. The simpletst one to 
compute is Mki|k 2 ^ 1 ~ §> • • • > ~ §) and the result is given as follows: 


-i 


^Ki|K ' 

Then the Sklyanin measure can be determined as (Ojk = Oj — 9k) 

n j<k(yj - yk) 


MKi|K 2 ' 


X = 


(& 12 ) L Y[{8jk)(8jk + i)(8jk ~ i) Res 

j<k 


y=x 


-ITn Qe ( yn)Qg (ijr 


(A.6) 


(A.7) 


In fact, this measure factor correctly reproduces the result given in [50], in which the overall 
normalization factor was determined by comparison with the known formulas. 


B The vacuum projection 


In this appendix we prove the formula (2.37) 


/ki|K 2 ( x ) 


(«i r ~( d 2 ) 


N* 


(B.l) 


where /ki|k 2 ( x ) = Ki|k 2 ( x I t L ) is the projection of the SoV basis on the vacuum state. 
The numbers N± are the numbers of pluses and minuses in the state k 1 |k 2 ( x I an d they are 
given in (2.22). We start with the explicit expression of /k 1 |k 2 ( x )) 


/ki|k 2 ( x ) 


l 


(t L in 

?c=i 


RlCiK 2 |o(^fc ) 

QoK) 



(B.2) 


It will be useful to write 


gK 2 1 |t i ) = ^e-“ 2S_ |t L ), (t L | = (t L |e a2S ", a 2 = c 2 /d 2 . (B.3) 

When all Sk = +, the expression (B.2) becomes 

/ Kl |K 2 (+,...,+) = ( t L I Skj I | L ) = ( d 2 ) L . (B.4) 

When the raising operators A Kl k 2 |o(^) are present we are going to insert the identity 
1 = e~° 2 s e + “ 2 s between consecutive operators, then we compute 

e a2S ~A KlK2l0 (8+)e~ a2S ~ = e a * s ~ (a 12 A* + b 12 C*) e~ a2S ~ . (B.5) 

Here and in the following we use for simplicity the notation A J = A{9j +i/ 2), etc. For this 
purpose we use the formula 

e a0 102 e -aOi =c>2 + a[0l ,0 2 ] + ±a 2 lO u [O u 0 2 ]] + ..., (B.6) 

and the commutators 

[S~Mu)] = D(u) - A(u) , [S~,A(u)] = B(u) , [S~,D(u)] = -B(u ), (B.7) 
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so that 


^KiK 2 |o(^) = e0l2S ^KiK 2 |o(^ i") e ° 2S ( B - 8 ) 

= 012 {A j + a-2 Bi) + b 12 ( C j + a 2 (D> - A’) + o| R J ) . 

In the next step we use that D 7 It 1 ') = C - 7 | | L ) =0 and A J | t L ) = QeCfy + '<) | t L )• 
Therefore, from the rightmost factor we obtain 

Q +(«+) 5 5 11 1 ) = ((«12 - 0 ^ 12 ) - (01202 + 6,2 “ 2 )Q #( ^ + i)) ' ^ ’ ' (R9) 

The first term is what we need to obtain (B.l) since 

ai 2 - a 2 b\ 2 = aia 2 + bic 2 - ^■(a 1 b 2 + bid 2 ) = ^ . (B.10) 

d 2 d 2 

The unwanted second term in (B.9) has to be commuted with the A, B 1 C, D terms from 
the next factors, and finally act on (\ L | on which it will vanish. The commutators will 
mostly give terms which vanish on It 1 ')) as one can see from the algebra (the specific 
coefficients are irrelevant) 

C j B k ~ B k C j + A k D j - A j D k , (B.ll) 

I) 1 B k ~ B k D j + B j D k , 

A j B k ~ B k A j + B j A k . 

The only non-vanishing terms are those coming from the last line, and they reproduce B r 
s. Repeating the procedure recursively, and using that (t L \ B-> = 0, one gets the desired 
result. We also need the opposite overlaps /k^k^x) = (i L I x )ki|k 2 > expressed as 

L 

/ki 2 ( x )= u L igK 2 n 

k= 1 

The last equality can be proven as above. Again, we write 

| l L ) = e - Q 2 s | 4 . l ) ; and (i L | g Ka = (d 2 ) L (\, L | e ° 2 s . (B.13) 

First we consider the case all sk = —, which gives 

ik 12 (- ..., -) = (| L I gK 2 || L ) = ( d 2 ) L . (B.14) 

When the creation operators are present we use again 

A Ki2]o ( 0J) = e a * s ~A Kl K 2 | O (0-)e-“ 2S " = (B.15) 

= 012 (Aj + ol 2 Bj ) + 612 (Cj + 02 (Dj — Aj ) T- o| Bj ) . 


^KiK 2 |o(^ 
Qe (®k ) 


11^ ) = ( a i) N + (d 2 ) 


\N X 


(B.12) 
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where now Aj = A(9J), etc. Furthermore, since Dj \ \. L ) = Bj \ \. L ) = 0 and Aj \ \. L ) = 
Qe(Qj) I | L )> the unwanted terms contain Cj at the right and they can be shown to vanish 
by using the commutation relations 

BjCk ~ CkBj + A.jDk — AkDj , (B.16) 

DjCk ~ CkDj + CjDk , 

AjCk ~ CkAj + CjAk . 

Finally, in order to compute the overlaps between the SoV bases of a whole chain and two 
subchains we need the initial condition 

<K X ) = Ki|o(+) ' ' ' > +| ® 0|K 2 ( _ > ' ' ' > “I I x )k!|K 2 = (t Li I ® (t Li I l x )Ki|K 2 = (t L l x )Ki|K 2 

_!A.L\„ TT ^KiK 2 |o(^) 2 , ,L\ ,0 1-71 

_<t lBt *5r«jwn u> ' <B17) 

Now we use 

|) = e _/?2 5 || L ), and (f L | g Ka = (b 2 ) L (i L \ e p2 s /3 2 = a 2 /b 2 . (B.18) 

Playing the same game as before we obtain that the contribution of each raising operator 
amounts to 

s j 4 Ki2 | 0 (Qj ) e ^ 2 s = (ai 2 — ^ 2 ^ 12 ) Aj + ... = ——Aj + ... , (B.19) 

so that 

<Kx) = (~h) N + (b 2 ) N - . (B. 20 ) 

It is reassuring to see that none of the various coefficients related to the SoV basis depend 
on ci and C 2 , which allows one to work with upper triangular twist matrices. 

C The splitting function 

In this subsection, we derive the difference equation for <h(y;y 2 |x) and give a solution in 
terms of T-functions,. Let us consider the following quantity 

(yi;y2|-B Kl |K 2 (w)l x >, (C.i) 

where we have omitted the indices indicating the twist of the spin chains. We first act the 
operator B Kl \ K2 (u) on the right, since |x) is the eigenstate of the double twisted F>-operator, 
we have 

L 

(yi;y2| J BK 1 |K 2 («)l x > = &12 - ^)(yi;y2l x >- (C.2) 

fc=i 
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On the other hand, we can act the operator B Ki |k 2 ( u ) on th e left by using the following 
relation 


- B Ki|K 2 ( u ) — ^Ki|o( w )- B 0|K 2 ( u ) + -®Ki|o( m )A)|K 2 ( u )- (C-3) 

By taking u = y\^ and u = y 2 ,j, we can derive two sets of difference equations 

&i2Qx(yij)$(yi;y2|x) = b 2 Q y2 (yij)Qg 1 (yij)^(-■ ■ ,yij -*,■■■ ;y2|x) (C. 4 ) 

bi 2 Q x (y 2 ,j)^(yi;y 2 |x) = b ] Q y] (y 2 j) Qe 2 (y 2 j)^(yi; ■ • • ,y 2 j + v • • l x ). 

Before solving these equations, let us note that the above two sets of difference equations 
(C.4) do not shift Xk- This implies that when we solve the above equations, we treat Xk 
as constant and the solution of these equations can only fix <3?(yi;y2|x) up to an initial 
condition <^(x) defined in equation (B.17) and computed in appendix B 

<Kx) = {~bi) N + (b2) N ~ ■ (C. 5 ) 

The strategy of solving equations of type (C.4) is to assume that the solution takes fac¬ 
torized form and decompose the equations into simpler difference equations which can be 
solved straightforwardly in terms of T-functions. Let us start by assuming that the solution 
of the difference equation takes the following form 

0 s(yi;y 2 |x) = p(yi;y 2 |x)^(yi;y 2 |x) (C. 6 ) 

where /o(yi;y 2 |x) satisfies the equations 

b 2 

p( yi;y2|x) = — Q y2 (yij)p(- ■ ■ ,yij ;y2|x) (C. 7 ) 

012 

p(yi;y2|x) = p- Q yi (y2j) p(yr, ■ ■ ■ ,y2j + *,••• |x). 

012 

The function y 2 |x) has to take care of the remaining factors in the difference equations 
(C.4) and satisfies the following equations 

Qx(?/ij)¥>(yi;y2|x) = • • ,yij ;y 2 |x) (C. 8 ) 

QAv2,j) 7 , (yi;y2|x) = Qg 2 (y2,j)<p{yv, ■ ■ ■ ,y2j + i, ■ ■ • |x). 

It turns out the solution of equations in (C.8) has the following factorized form 

Li l 2 

<7 , (yi;y2|x) = JJ^i(yi,jl x ) JJ <^2 (2/2 j|x) (c. 9 ) 

3 = 1 i =1 

where (pi(y\x) and <£> 2 (, 2 |x) satisfy the following equations 

Qx(y)^i(y|x) = Qg^yfyiiy - *|x), (C.10) 

Q x {y) P2{z |x) = Qg 2 (z)tp2{y + *|x). 
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The solution is in the end given by (2.37): 3>(yi;y2|x) = twist x Gamma, with 


twist = (ib 12 /h) N + (~ib 12 /b 2 ) N - (-h) N + (b 2 ) N ~ , 

_ r (i{0+ - 02)) r (1 - z(yi - flr)) £ (1 + i(x - 0 +)) r (1 - i(x - 07 2 )) 
r (*(yi - y 2 )) r (1 - *(y 2 - 0 2 )) r (1 + i(x - yi)) T (1 - i(x - y 2 )) ' 
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